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ABSTRACT 


In  this  study,  the  author  develops  a  MATLAB  simulation  of 
area  search  with  acoustic  sensors  modeled  by  the  Poisson  Scan 
model  and  the  Lambda-Sigma  (/I— cr)  model.  Detection  time  results 
are  compared  to  those  given  by  the  much  simpler  Random  Search 
formula.  Random  Search  was  found  to  closely  approximate  the  more 
complex  models  if  detection  range  was  selected  correctly. 
Guidelines  for  selecting  the  Random  Search  detection  range  were 
developed . 
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EXECUTIVE  SUMMARY 


In  this  study,  the  author  develops  a  MATLAB  simulation  of 
area  search  with  acoustic  sensors  modeled  by  the  Poisson  Scan 
and  the  Lambda-Sigma  (/I  —  cr)  models. 

Both  the  Poisson  Scan  model  and  Lambda-Sigma  (1  —  cr)  model 
simulation  results  are  found  to  be  approximately  exponentially 
distributed,  which  is  consistent  with  the  Random  Search  model. 
Thus,  Random  Search  with  the  proper  deterministic  detection 
range  R  can  be  used  to  closely  approximate  simulation  results 
obtained  using  the  Poisson  Scan  and  Lambda-Sigma  (1  —  cr)  models. 

Both  the  Poisson  Scan  and  Lambda-Sigma  (1  —  cr)  detection 
time  results  vary  with  A,  cr,  and  i?(50).  The  author  constructs  a 
regression  model  based  on  the  simulation  results  and  finds  an 
approximate  linear  relationship  between  i?(50)  and  the  best-fit 
R  for  reasonable  A  and  cr.  The  gradient  of  the  regression  line 
depends  on  the  values  of  A  and  cr  .  Thus,  it  is  possible  to 
estimate  the  best-fit  R  from  problem  parameters  A,  cr,  and  i?(50) 
and  then  use  Random  Search,  instead  of  Monte  Carlo  simulation, 
to  predict  the  effectiveness  of  area  search. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Random  Search  is  one  of  the  most  well-known  and  used  models 
for  area  search.  However,  because  of  its  simplicity,  it  also 
has  inherent  limitations.  For  example,  it  assumes  that  the 
searcher  and  target  search  areas  are  identical  and  that  the 
searcher  uses  a  perfect  cookie-cutter  sensor.  Captain  Gi  Young 
Kim  1  ,  ROK  Air  Force,  studied  these  issues  previously  and 
suggested  generalizations  to  the  Random  Search  formula  for 
cases  where  the  searcher  and  target  areas  are  not  identical. 
However,  there  are  remaining  issues  to  address,  particularly 
how  well  the  Random  Search  model  can  approximate  probabilistic 
sensors,  such  as  those  modeled  by  the  Poisson  Scan  and 
Lambda-Sigma  (A  — a)  models. 

B.  RESEARCH  QUESTIONS 

The  primary  research  questions  are  the  following: 

1.  When  using  the  Poisson  Scan  model  and  the  Lambda-Sigma 
(A  — a)  model,  will  detection  times  be  approximately 
exponentially  distributed,  as  predicted  by  the  Random 
Search  model? 

2.  Can  a  non-cookie  cutter  sensor  be  replaced  with  a 
cookie-cutter  sensor  with  the  same  sweep  width  and 
maintain  approximately  the  same  search  performance? 

3.  Are  there  significant  differences  in  the  detection  results 
generated  by  simulations  using  the  Poisson  Scan  and 
Lambda-Sigma  (A  — a)  detection  models? 


^  Gi  Young  Kim,  "Development  and  Testing  of  a  New  Area  Search  Model  with 
Partially  Overlapping  Target  and  Searcher  Patrol  Areas,"  Master's  thesis. 
Naval  Postgraduate  School,  Monterey,  California,  2008. 
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C.  THESIS  ORGANIZATION 

The  thesis  is  organized  as  follows:  Chapter  II  reviews  the 
Poisson  Scan  and  Lambda-Sigma  (/I  —  cr)  acoustic  detection  models; 
Chapters  III  and  IV  present  a  MATLAB  simulations  of  these  two 
models;  Chapter  V  develops  the  best-fit  R  ,  which  is  the 
equivalent  cookie-cutter  detection  range  for  the  Poisson  Scan 
model  and  Lambda-Sigma  (1  —  cr)  model;  Chapter  VI  summarizes  all 
results  and  recommends  future  studies. 
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II.  OVERVIEW  OF  ACOUSTIC  DETECTION  MODELS 


A.  BACKGROUND 

In  this  chapter,  we  review  the  Passive  Sonar  Equation (PSE) 
and  use  it  to  develop  two  models  for  passive  acoustic  detection 
—  the  Poisson  Scan  model  and  the  Lambda-Sigma  (A  — a)  model. 

B.  THE  PASSIVE  SONAR  EQUATION 

The  passive  sonar  equation^  is  an  audit  of  energy  flow 
between  a  source  and  a  sonar  receiver.  All  units  are  in 
decibels (dB)  . 

The  Passive  Sonar  Equation (PSE)  is 

DNL  -DT  =  SE, 

where 

Ls  —  signal  level  of  source, 

propagation  loss  between  the  source  and 
receiver, 

DNL  =  detection  noise  level  at  the  receiver, 

DT  =  detection  threshold  =  signal  level  at  receiver 
required  to  achieve  a  of  0.5, 

SE  —  signal  excess. 

The  left  side  of  Passive  Sonar  Equation  is  the  signal 
available  at  the  receiver,  and  right  side  is  the  signal  required 
for  a  probability  of  detection (P^ )  of  0.5. 


2  The  references  for  this  section  are  the  unpublished  lecture  notes  of 
Professor  James  N.  Eagle,  "Acoustic  Detection  Models,"  2009. 


3 


1. 


Figure  of  Merit  (FOM)  and  i?(50) 


Figure  of  merit  (FOM)  and  F(50)  are  used  as  measures  of 
ef f ectiveness  (MOF)  for  sonar  detection  performance.  FOM  is 
the  maximum  propagation  loss  a  sonar  system  can  absorb  in  dB  and 
still  produce  a  0.5  probability  of  detection.  i?(50)  is  the 
maximum  range  between  the  source  and  receiver  resulting  in  a 
0.5  probability  of  detection. 

FOM  =  Ls-  DNL  -  DT 

=  the  transmission  loss  between  source  and  receiver 
resulting  in  P(i=0.5. 

R(S0)  =  the  maximum  range  where  Ls  =  FOM . 

2.  Acoustic  Modeling  Assvimptions 

•  Mean  SE  =  SE  —  FOM  —  .  That  is,  the  Passive  Sonar 

Equation  is  used  to  calculate  mean  signal  excess  at  the 
receiver . 

•  SE~N(SE,a^),  where  cr  is  typically  3-9dB .  We  assume  the 
signal  excess  is  a  normally  distributed  random  variable 
with  the  mean  given  by  the  Passive  Sonar  Equation,  and 
a  is  specified. 

•  Detection  occurs  if  and  only  if  SE  >  OdB  .  This  is  a 
threshold  crossing  model.  Detection  occurs  when  the 
random  SE  is  nonnegative  and  at  no  other  times.  And 
probability  of  detection  (called  here  instantaneous 
probability  detection,  IPD)  is  P(SE>0). 
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Figure  1 .  IPD  when  SE  < 0 


3 .  Implications 

1PD(SE)  =  PiSE  >  0),  where  SE~N(SE,a^) 

=  0(SE/a)  . 

where  0  is  the  cumulative  normal  distribution  function.  This 
relationship  between  SE  and  the  probability  of  detection  is 
illustrated  in  Figure  1. 

C.  THE  FIXED  SCAN  AND  POISSON  SCAN  MODELS 
1 .  Fixed  Scan  Model 

The  Fixed  Scan  model ^  assumes  that  probabilistically 
independent  detection  opportunities  occur  at  regularly  spaced 
time  intervals.  It  is  attractive  analytically,  since  there  is 
a  simple  formula  for  calculating  the  probability  of  detection 
over  any  time  interval  [0,  t] .  Let  S  be  the  set  of  detection 
opportunity  times  t  within  the  interval,  and  let 

/PD(t)  =  0(SE(T)/a). 


3  Alan  R.  Washburn,  Search  and  Detection,  4th  ed.  3-4p. 
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be  the  detection  probability  for  a  scan  containing  time  t,  where 
the  signal  excess  SE(t)  is  acknowledged  to  depend  on  time.  Then 
the  cumulative  detection  probability  over  the  interval  is 

P(T  <  t)  =  Frit)  =  1  -  Uresa  "  IPD(.t))  . 

In  other  words,  there  will  be  a  detection  at  [0,  t]  unless 
detection  fails  at  every  one  of  the  independent  scans  within 
the  interval.  However,  the  Fixed  Scan  model  suffers  from  its 
dependence  on  an  arbitrarily  selected  origin  of  time.  The 
Poisson  Scan  model  was  developed,  in  part,  to  address  this 
issue . 


2.  Poisson  Scan  Assvimptions 

•  Independent  detection  opportunities  occur  at  Poisson 
times  with  rate  1 (units  :  1/time) . 

•  Probability  of  detection  at  time  t  ,  given  that  a 
detection  opportunity  exists  at  t,  is  I  PD  it)  =  0iSEit)/a). 

3 .  Development 

From  the  Poisson  Scan  assumption,  the  instantaneous 
detection  rate  at  time  t  is 


Yit)  =  A0iSEit)/a). 

Then  it  follows  from  the  properties  of  the  non-homogeneous 
Poisson  Process  that 

PiT  <  t)  =  Frit)  =  1  -  exp  i-A  fl^0iSEit)/a)ds). 
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In  this  model,  1/1  is  the  mean  time  between  independent 
detection  opportunities  (glimpses) .  Probabilistic  independence 
is  assumed  for  each  glimpse,  with  1  specifying  the  average 
glimpse  rate.  By  varying  1,  the  sonar  can  be  made  as  effective 
(1  larger)  or  ineffective  (1  smaller)  as  desired,  subject  to 
sensor  or  operational  constraints.  The  sonar  system  can  only 
call  detection  at  the  discrete  glimpse  times,  and  otherwise  is 
assumed  to  be  processing  previously  received  data. 

D.  THE  LAMBDA- SIGMA  (A -a)  MODEL 

In  the  Poisson  Scan  model,  detection  can  occur  at  time 
t  only  if  a  detection  opportunity  occurs  at  t,  and  an  independent 
draw  from  SE  ~  N(^SE,  a^)  is  greater  than  or  equal  to  0 . 

A  potentially  more  realistic  model  would  assume  that  5£'(t) 
is  a  random,  continuous  function  of  time  t  (that  is,  a 
continuous  -  time  stochastic  process) ,  and  that  detection 
occurs  at  any  time  t  where  5£'(t)>0. 

1.  (1  — <7)  Model  Assiamptions  and  Development 

Let 

SE(t)  ^EE(t) 


where  SE(t)  is  the  mean  signal  excess  at  time  t  ,  which  is 
computed (as  in  the  Poisson  Scan  model)  from  the  PSE  as 


SE(t)  =  EOM(t)  -iVi^(t). 

X(t)  is  a  0-mean  stochastic  process.  More  specifically,  X(t)  is 
a  (A  — a)  jump  process,  described  below. 
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Each  random  sample  path  of  X(t)  is  a  step  function,  where 
the  duration  of  each  step  (the  time  between  jumps)  is 
exponentially  distributed  with  mean  1/A,  and  the  height  of  each 
step  (the  value  of  X(t))  is  normally  distributed  with  mean  0  and 
variance  cr^  . 


Figure  2.  (A  — a)  Detection  Process  Example 

As  illustrated  in  Figure  2,  detection  occurs  whenever 
SE(t)  >  X(t)  .  X(t)  can  be  thought  of  as  the  time-varying  sum  of  all 
random  components  of  FOM  and  . 

2.  Computing  when  SE  is  Constant 

Assume  SE(^t)  =  SE  —  k,  t  E  [0,  t]  .  We  will  develop  an  expression 
for  the  probability  of  initial  detection  occurring  between  time 
0  and  t,  which  we  call  Fj{E) . 
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Frit)  =  P(X(t)  <  k  ,  for  some  t  G  [0,  t]) 

=  1  —  P{X{t)  >  k,  for  all  t  G  [0,  t]) 

=  1  —  T,i=oP(,X(t)  >  k,  for  all  t  G  [0,  T]|i  jumps  occur  in  [0,  t]) 
X  P(i  jumps  occur  in  [0,  t]) 


=  1  -  Sr.o(l  -  4- 


=  1  -  [1  -  4.(k/<r)]e-->tsr.0 


ai-0{k/a)]Aty 

a 


=  1  -  [1  -  since  =  e 

=  1  -  [1  - 

=  1  -  [1  -  . 


A 


r(<k'tc(tioii  at  tiiiio  0)  =  <IHSE/(t} 

y _ 


time  t 


Figure  3.  (X  —  a)  PDF  and  CDF  Functions  for  Time  of 

Initial  Detection  T  when  SE  is  Constant 
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At  time  0,  the  probability  of  detection  is  0(SE/a) . 

The  mean  time  to  detection,  E(T),  and  the  density  function, 
for  time  of  initial  detection  are 


E(T)  = 


1  -  0iSE/a) 
A0(SE/a)  ' 


frit)  =  X0iSE/o)[l  -  for  t  >  0. 


Also,  Pit  —  0)—  0iSE/a). 


Note  that  this  probability  distribution  has  a  discrete  portion 
(r  =  o)  with  units  of  probability  (unitless)  and  a  continuous 
portion  (T  >  0)  with  units  of  probability/ time  (or,  1/time). 
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III.  SIMULATION  OF  POISSON  SCAN  MODEL 


A.  DESCRIPTION  OF  POISSON  SCAN  SIMULATION  MODEL 

1.  Characteristics  of  the  Searcher 

The  searcher  is  not  allowed  to  search  for  a  target  outside 
the  search  area.  In  addition,  the  searcher's  initial  position 
is  uniformly  distributed  inside  of  the  search  area.  After  that, 
the  searcher  selects  his  course  randomly,  independent  of  the 
target's  movement.  The  course  change  event  is  determined  by 
Poisson  process  with  rate  I5 .  Independent  scan  opportunities 
also  occur  according  to  a  Poisson  process,  but  with  rate  .  Also, 
it  is  assumed  that  signal  transmission  loss'^  follows  either  the 
spherical  spreading  law  (iV^  =  20  logio  r)  or  the  cylindrical 
spreading  law  (iV,^  =  10  log^o  •  The  searcher  speed  was  fixed  at 
15  nm/hour . 

2.  Characteristics  of  the  Target 

Like  the  searcher,  the  target  is  not  allowed  to  move 
outside  of  the  search  area  and  has  an  initial  position  uniformly 
distributed  over  the  search  area.  The  logic  of  the  target 
movement  is  the  same  as  that  of  the  searcher;  that  is,  the  target 
has  its  own,  independent  Poisson  process  with  course  change  rate 
If.  The  target  speed  was  fixed  at  5nm/hour. 
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Robert  J.  Urick,  Principles  of  Underwater  Sound,  3rd  ed.  lOlp. 
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B.  COMPUTER  ALGORITHM 

1.  Input  [Units] 

•  Number  of  simulation  replications,  N^gpg  =  SOO . 

•  Maximum  simulation  time,  t^ax  =  lS0[hour] . 

•  The  length  of  search  area  in  X  direction,  =  150[nm] . 

•  The  length  of  search  area  in  Y  direction,  ly  =  lS0[nm] . 

•  Searcher  speed,  V  [nm/hour] . 

•  Target  speed,  U  [nm/hour], 

•  Searcher's  scan  rate,  Ai  [Glimpses /hour] . 

•  Searcher's  course  change  rate,  Ag  [1/hour] . 

•  Target's  course  change  rate.  If  [1/hour] . 

•  The  unit  time  of  simulation.  At  [hour] . 

•  The  size  of  search  area,  Ag  =  l^-  x  ly  . 

•  Figure  of  Merit,  FOM  [dB] . 

•  Transmission  loss,  =  20  \og(^Distance)  [dB] . 

•  Signal  Excess,  SE  =  FOM  —  N^^. 

•  Variance  of  Signal  Excess,  cr  [dB] . 

2.  Functioning  of  the  Poisson  Scan  Simulation 

When  a  new  replication  begins,  the  initial  positions  of 
the  searcher  and  the  target  are  chosen  from  a  uniform 
distribution  over  the  search  area.  The  initial  headings  are  also 
drawn  uniformly  between  0  and  2n  radians.  The  subsequent  course 
changes  for  the  searcher  and  target  occur  according  to  Poisson 
processes  with  rates  Ag  and  Af-  ,  respectively.  According  to 
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Captain  Gi  Young  Kim' s  thesis^,  the  recommended  value  for  the 
course  change  rate  is  A  =  2V/^  ly  ,  which  implies  that  on 
average,  two  course  change  events  occur  during  the  time  required 
for  the  searcher  to  go  from  edge  to  edge  in  the  search  area. 

When  the  searcher  or  target  encounters  an  area  boundary, 
a  random  reflection  occurs.  After  a  reflection,  the  new  course 
in  radians  is  Uniform_Random(^Cy  —  O.S,Cy  +  O.S),  where  Cy  is  the 
perpendicular  course  from  the  reflection  boundary.  This  scheme 
was  recommended  and  tested  by  Captain  Gi  Young  Kim®  to  achieve 
an  approximate  uniform  distribution  of  the  searcher  and  target 
tracks  in  the  search  area. 

In  the  Poisson  Scan  model,  there  are  two  ways  to  simulate 
detection  times.  One  uses  the  instantaneous  probability 
detection  at  specific  time  t,  which  follows  the  Poisson  process 
with  rate  1.  This  is  because  in  the  Poisson  Scan  model,  detection 
opportunities  occur  according  to  a  Poisson  process  with  rate 
1.  So,  the  mean  time  between  independent  detection  opportunities 
is  1/1.  The  distance  between  the  searcher  and  the  target  is 
calculated  at  specific  time  t.  Then,  the  mean  signal  excess  (SE) 
is  determined  by  using  the  Passive  Sonar  Equation (PSE)  allowing 
the  calculation  of  IPD(SE)  .  The  next  step  is  generating  a 
standard  uniform  random  number  and  comparing  it  with  the  IPDiSE) 
If  Uniform  Random{Q,E)  <  IPD(^SE) ,  then  a  detection  occurs  at  that 
specific  time  t. 


^  Gi  Young  Kim,  "Development  and  Testing  of  a  New  Area  Search  Model  with 
Partially  Overlapping  Target  and  Searcher  Patrol  Areas,"  Master's  thesis. 
Naval  Postgraduate  School,  Monterey,  California,  2008. 

6  Gi  Young  Kim,  KOREA  AIR  FORCE  CAPTAIN,  graduated  from  the  Naval 
Postgraduate  School  in  2008. 
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Figure  4.  Simulation  when  calculating  distance  and  SE  at 

specific  time  t 


The  other  way  to  simulate  detections  is  to  check  for  a 
detection  at  each  At.  The  distance  between  the  searcher  and  the 
target  is  calculated  at  every  time  step  At,  and  the  mean  signal 
excess  is  computed. 


After  that,  we  compute  the  instantaneous  probability  of 
detection  IPD(SE)  —  <P(SE(t)/a)  and  detection  rate  /(t)  =  A<?(5£'(t)/(T)  . 
At  each  time  step  At,  if  Uniform  Random(0,l)  <  ,  then  a 

detection  occurs. 
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Figure  6.  Frit)  for  different  simulations 

As  illustrated  in  the  simulation  results  of  Figure  6,  the 
results  of  these  two  simulation  methods  are  very  similar,  and 
theoretically  should  be  identical. 

For  this  thesis  we  used  the  At  simulation  method.  In  order 
to  approximate  a  continuous  simulation.  At  should  be  small. 
However,  too  small  a  At  results  in  too  many  calculations.  In  the 
next  section,  we  answer  the  question,  "How  small  should  At  be 
to  produce  accurate  results  without  creating  excessively  long 
simulation  runs?" 
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Figure  7 .  Frit)  for  various  values  of  At 


3 .  Output 

Extensive  simulation  experimentation  showed  that  At  — 
was  the  largest  At  resulting  in  reliable  and  repeatable  results 
Thus  we  move  the  searcher  half  the  distance  of  the  i?(50)  at  each 
time  step.  Typical  results  are  shown  in  Figure  7. 

In  order  to  find  best-fit  R,  which  is  the  equivalent  range 
of  a  cookie-cutter  sensor  to  the  Poisson  Scan  Model,  the  author 
experimented  with  various  values  of  1,  cr,  and  FOM . 
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FOM  =  90dB 


FOM  =  87dB 


R  =  18.15n?n  y  =  0-0250 
FOM  =  8AdB 


R  =  13.08n?n  y  =  0-0181 
FOM  =  80dB 


0.9 


0.8 

0.7 

0,6 

Q. 

Q 

0.5 

U 

0,4 

0,3 

0.2 

0.1 


0 


Time 


R  =  8.24nm  y  =  0.0114 
FOM  =  77  dB 


FOM  =  74:dB 


Figure  8 . 


Frit)  for  various  FOM  in  Poisson  Scan  model 
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Figure  9.  Frit)  for  various  Sigma  in  Poisson  Scan  model 
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igure  10.  ^rCO  various  Lambda  in  Poisson  Scan  model 
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The  Green  line  represents  the  Poisson  Scan  model  formula. 


P(T  <  t)  =  =  1  —  exp  (—1  I  0(SE(s)/a)ds). 

Js=o 

The  red  line  represents  the  formula  resulting  from  a  random 
search  with  a  cookie-cutter  sensor  that  has  best-fit  R  range. 

P(T  <  t)  =  Frit)  =  1  -  exp  (-2RVt/As), 


5  _  —  ln(l  —  Fjjty^ax))  * 

~  2Vt 

^  ^  ^max 

where  V  is  the  approximate  mean  relative  speed  between  target 
and  searcher  determined  as  follows: 

V^-\  m  -2VUcos9d9. 

TTJo 

The  blue  line  represents  the  results  of  the  simulation. 

The  results  show  that  for  the  model  parameters  examined 
here,  the  Poisson  Scan  model  gives  detection  times  that  are 
approximately  exponentially  distributed,  as  does  the  Random 
Search  model.  The  conclusion  is  that,  the  Random  Search  model 


using  the  proper  detection  range  R  can  closely  approximate 
detection  results  given  by  the  Poisson  Scan  model. 

A  potentially  more  accurate  procedure  for  determining  R 
would  be  to  plot  — In  (1  —  F7’(t))  versus  t  and  to  solve  for  the 
best-fit  slope  y.  Then, 
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IV.  SIMULATION  OF  LAMBDA-SIGMA  (A  -  fl)  MODEL 


A.  DESCRIPTION  OF  LAMBDA-SIGMA  (A  -  ff)  SIMULATION  MODEL 

The  Lambda-Sigma  (A  —  cr)  detection  model^  is  distinguished 
from  the  Poisson  Scan  model  by  allowing  detection  at  any  time 
during  the  scan,  rather  than  only  at  specific  scan  times.  This 
model  assume  that  SE(t)  is  a  random,  continuous  function  of  time 
t.  Then, 


SEit)  =  SEit)  -  A(t), 

where  X(t)  is  a  0-mean  stochastic  process  where  the  duration  of 
each  step  is  exponentially  distributed  with  mean  1/A,  and  the 
height  of  each  step  is  normally  distributed  with  mean  0  and 
variance  cr^ .  Detection  occurs  at  any  time  t  where  5£'(t)>0. 


The  references  for  this  section  are  the  unpublished  lecture  notes  of 
Professor  James  N.  Eagle,  "Acoustic  Detection  Models,"  2009. 
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B.  COMPUTER  ALGORITHM 

1.  Input  [Units] 

•  Number  of  simulation  replications,  N^gpg  =  SOO . 

•  Maximum  simulation  time,  t^ax  =  lS0[hour] . 

•  The  length  of  search  area  in  X  direction,  ly.  =  150[nm] . 

•  The  length  of  search  area  in  Y  direction,  ly  =  150[nm] . 

•  Searcher  speed,  V  —  IS  \nm/hour^ . 

•  Target  speed,  U  =  S  [nm/hour] . 

•  Searcher' s  course  change  rate,  Ag  =  AV/yJly  X  ly  [1/hour] . 

•  Target's  course  change  rate,  —  2U / ^ ly  X  ly  [1/hour]. 

•  The  unit  time  of  simulation.  At  =  i?(50)/2  X  K  [hour]. 

•  The  size  of  Search  area,  Ag  =  ly  x  ly  . 

•  Figure  of  Merit,  FOM  [dB] . 

•  Transmission  loss,  =  20  \og(^Distance)  [dB] . 

•  Signal  Excess,  SE  =  FOM  —  N^^. 

•  Jump  rate  of  stochastic  process,  Xi  [times/hour]. 

•  Standard  deviation  of  jump  in  stochastic  process,  a[dB] . 

2 .  Functioning  of  the  Program 

With  the  Lambda-Sigma  (1  —  cr)  model,  the  duration  of  each 
step  is  exponentially  distributed  with  mean  lA,  and  the  height 
of  each  step  is  normally  distributed  with  mean  0  and  variance 
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Figure  11.  Example  of  (A -a)  model  simulation 

In  Figure  11,  the  blue  line  represents  one  realization  of  the 
0-mean  stochastic  process  X(t),  and  the  red  line  represents  mean 
signal  excess.  At  each  time  step  At,  if  (t)  >  A'(t)  ,  then  a 
detection  occurs.  In  the  Figure  11,  the  first  detection  occurs 
at  55  hours. 

3 .  Output 

In  order  to  find  the  best-fit  R,  which  is  the  equivalent 
range  of  a  cookie-cutter  sensor  to  that  of  the  Lambda-Sigma 
(A -a)  model,  the  author  experimented  with  various  values  of 
model  parameters  A,  cr,  and  FOM . 
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Figure  12 .  Frit)  for  various  FOM  in  (1  —  cr)  model 
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Figure  13. 


Frit)  for  various  Sigma  in  (1  —  cr)  model 


R  =  S.03nm  y  =  0-0069 


R  =  S.SSnm  y  =  0-0076 


Lambda  =  l.S[l/hour] 


Lambda  =  2[l/ftour] 


R  =  S.86nm  y  =  0.0079 


Lambda  =  2.5[l/ftour] 
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R  =  5.89nm  y  =  0.0081 


Lambda  =  3[l/ftour] 


R  =  5.99nm  y  =  0.0082 


R  =  6.06nm  y  =  0.0083 


Figure  14 . 


F^Ct)  for  various  Lambda  in  (1  —  cr)  model 
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The  red  lines  in  Figures  12-14  show  Random  Search  results 
with  a  cookie-cutter  sensor  having  best-fit  R  range.  The  blue 
line  represents  the  results  of  the  simulation.  As  was  the  case 
with  the  Poisson  Scan  model,  simulation  results  using  the 
Lambda-Sigma  (1  —  a)  model  showed  an  exponential  time  to  initial 
detection.  Thus,  Random  Search  with  the  correct  deterministic 
detection  range  R  can  be  used  to  closely  approximate  simulation 
results  obtained  using  the  stochastic  Lambda-Sigma  (A  —  a)  model 
The  task  that  remains  is  to  estimate  from  problem  parameters 
the  appropriate  value  of  R. 
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V.  ESTIMATION  OF  BEST-FIT  R 


A.  LATERAL  RANGE  CURVE 

1 .  Description  of  the  Lateral  Range  Curve 


In  order  to  estimate  the  best-fit  R,  the  author  simulated 
lateral  range  curves  for  sensors  with  the  same  parameters  as 
the  Poisson  Scan  model  and  Lambda-Sigma  (1  —  cr)  model.  Sweep 
width  W  is  defined  to  be  the  area  underneath  the  lateral  range 
curve . 

W  =  I  l{x)dx, 

‘'  —  00 

where  /(X)  is  the  lateral  range  curve.  In  the  cookie-cutter 
sensor  with  detection  range  R, 


l{x)  = 


therefore,  R  =  W/2 . 


l,x€  [-7?,  R] 
0,  otherwise 


r 


Figure  15.  Sweep  Width  Interpretation 
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2 .  Computer  Algorithm  of  Lateral  Range  Curve 

In  the  lateral  range  curve  simulation  geometry,  shown  in 
Figure  15,  target  lateral  ranges  are  uniformly  distributed  over 

-L/2  to  L/2  and  target  tracks  are  straight.  Except  for  this 
characteristic,  input  elements  in  each  lateral  range  curve 
simulation  are  the  same  for  both  the  Poisson  Scan  model  and 
the  Lambda-Sigma  (A -a)  model.  The  lateral  range  curve  is 
determined  by  the  ratio  of  the  number  of  detections  that  occur 
to  number  of  replications  at  each  lateral  range  from  the 
searcher . 


Figure  16.  Simulation  result  of  lateral  range  curve 

Then,  the  area  underneath  the  lateral  range  curve  is  sweep  width, 
which  can  be  computed  numerically. 
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3. 


Output  of  Lateral  Range  Curve 


The  following  tables  compare  the  half-sweep  width  computed 
from  lateral  range  curve  simulations  to  the  best-fit  R  obtained 
from  Poisson  Scan  and  Lambda-Sigma  (1  —  cr)  simulations,  for 
various  FOM ,  X,  and  cr  values  and  assuming  spherical  spreading 


FOM\dB] 

74 

77 

80 

84 

87 

90 

Poisson  Scan 
model [nm] 

1 .40 

2 .74 

4  .  94 

8 .24 

13.08 

18 . 15 

(1  —  0')  model[nm] 

2.80 

3.97 

5.37 

8.81 

12 .81 

19.07 

Lateral  range  half 
sweep  width  of  P-S 
model [nm] 

2 .14 

3 .75 

6.31 

10.26 

16.15 

28 . 42 

Lateral  range  half 
sweep  width  of 
(1  —  0')  model  [nm] 

3.07 

4.49 

6 . 65 

11.31 

16.87 

25.63 

Table  1.  Comparison  of  R  between  Poisson  Scan  model,  (A  -  o') 
model  and  half  sweep  width  of  each  model  when 

TJn  1  T  1  17  ir  _  nm 

O'  =  3dB,  A  =  2 ,V  =  15 ,U  =  5 

hour  hour  hour 


Sic]ma\dB] 

0 

1 

2 

3 

4 

5 

6 

7 

Poisson  Scan 
model  [nm] 

3.49 

3.51 

3.75 

4 . 68 

5.25 

6.52 

8.51 

11 . 83 

(1  -  (j)  model[nm] 

4.61 

5.01 

5.11 

5.89 

6.56 

7 . 84 

10.08 

12 . 68 

Lateral  range  half 
sweep  width  of  P-S 
model [nm] 

3.68 

4.29 

5.12 

6.31 

8.07 

10.62 

14.56 

20.78 

Lateral  range  half 
sweep  width  of 
(1  -  (j)  model  [nm] 

5.30 

5.55 

5.98 

6.64 

7.70 

9.36 

11.86 

15.42 

Table  2.  Comparison  of  R  between  Poisson  Scan  model,  (1  -  o') 
model  and  half  sweep  width  of  each  model  when 
FOM  =  80dB,  X  =  2—,V=  1S—,U  =  5  — 

hour  hour  hour 
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/l[l/ftour] 

0.5 

1 

1 . 5 

2 

2.5 

3 

Poisson  Scan 
model [nm] 

1.79 

2 .77 

3.86 

4.68 

4 . 94 

5.25 

(A  — a)  model[nm] 

5.03 

5.58 

5.86 

5.89 

5.99 

6.06 

Lateral  range  half 
sweep  width  of  P-S 
model [nm] 

0.89 

2 . 98 

4.89 

6.32 

7 .40 

8.23 

Lateral  range  half 
sweep  width  of 
(A  —  O')  model  [nm] 

5 .73 

6.07 

6.36 

6.64 

6.84 

7 .02 

Table  3.  Comparison  of  R  between  Poisson  Scan  model,  (1  -  o') 
model  and  half  sweep  width  of  each  model  when 
FOM  =  80dB,  a  =  3dB,  V  =  1S—,U  =  5  — 

hour  hour 

These  results  showed  that: 

•  The  best  fit  R  is  strongly  dependent  on  FOM ,  A,  and  cr  for 
both  the  Poisson  Scan  and  Lambda-Sigma  (1-a)  models . 

•  The  Lambda-Sigma  (1  —  cr)  best  fit  R  somewhat  exceeds  the 
Poisson  Scan  best  fit  R,  and  this  might  be  due  to  the 
Lambda-Sigma  (1  —  cr)  model  starting  with  a  positive 
probability  of  detection  at  time  0. 

•  The  best  fit  R  generally  increases  with  increases  in  FOM , 
X,  and  cr  for  both  the  Poisson  Scan  and  Lambda-Sigma  (1  —  cr) 
models 

•  The  lateral  range  curve  model  produces  half-sweep  width 
values  which  can  significantly  exceed  the  best  fit  R 
values  produces  by  the  Poisson  Scan  and  Lambda-Sigma  (1  —  cr) 
models 

It  is  not  clear  why  the  lateral  range  procedure  produced 
detection  range  estimates  not  matching  well  with  those  of  the 
Poisson  Scan  and  Lambda-Sigma  (1  —  cr)  models;  but  it  is  possible 
that  the  lateral  range  assumptions  (infinite,  straight-line 
paths  with  uniformly  distributed  closest  points  of  approach) 
were  not  well  enough  met  in  the  area  search  simulation. 
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B.  i?(50) 


1.  Description  of  i?(50) 

Another  way  to  estimate  the  best-fit  R  is  to  use  i?(50), 
which  is  the  maximum  range  between  the  searcher  and  target 
resulting  in  an  instantaneous  detection  probability  of  0.5.  In 
other  words,  i?(50)  is  the  maximum  range  where  mean  signal  excess 
SE  is  zero.  For  this  analysis,  the  author  assumes  that  signal 
transmission  loss  follows  a  mixture  of  spherical  spreading 
(iV^  =  20  logio  r)  and  cylindrical  spreading  (iV^  =  10  logio  r)  • 
Specifically,  we  use  =  ISlogigr.  The  author  also  assumes  that 
the  source  level  of  the  target  is  in  the  range  84  to  95dB , 
detection  noise  level  at  the  receiver  DNL  is  45  dB  ,  and 
detection  threshold  DT  is  -15dB.  Therefore,  by  the  passive 
sonar  equation,  FOM  is  in  the  range  54  to  65  dB . 

2 .  Compute  i?(50) 

With  the  author's  assumptions,  the  FOM  and  i?(50)  can  be 
calculated . 


SE  =  FOM  - 

=  FOM  —  15  logio  ^  • 


Setting  SE  =  OdB ,  =  0.5,  and  Inm  =  1852m,  we  determine  that 


i?(50)  = 


FOM 

lOT^ 

1852 


This  value  can  be  compared  to  the  best-fit  R  values  obtained 
from  the  Poisson  Scan  and  Lambda-Sigma  (1  —  cr)  simulations. 


33 


3. 


Output  of  i?(50) 


For  acoustic  detection  models,  a  reasonable  1  is  1  to  2hr~^ 
and  a  reasonable  cr  is  2  to  AdB .  The  author  used  regression  to 
examine  how  well  the  computed  R(S0)  values  estimated  the 
best-fit  R  values  obtained  from  the  Poisson  Scan  and 
Lambda- Sigma  (1-a)  simulations . 


Figure  17 . 


Regression  result  of  Poisson 


Scan  model 
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Figure  18.  Regression  result  of  Lambda-Sigma  (1-a)  model 

As  shown  in  Figures  17  and  18,  the  relationship  between 
i?(50)  and  best-fit  R  is  approximately  linear.  In  addition,  the 
slope  of  each  regression  line  strongly  depends  on  the  values 
of  1,  cr,  and  the  model  being  used. 


1  [/ir-i] 

2  [/ir-i] 

2dB 

R  =  0.85  xi?(50)  -  1.71 

R  =  1.10  X  R(SQ)  -  1.67 

3dB 

R  =  1.06  xi?(50)  -2.17 

R  =  1.44  X  i?(50)  -  2.41 

AdB 

R  =  1.49  x«(50)  -3.05 

R  =  2.13  x«(50)  -3.79 

Table  4.  Regression  result  of  Poisson  Scan  model  at  each  A  and 

a 
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2[/ir-i] 

2dB 

R  =  1.05  xi?(50)  -0.16 

R  =  1.14  X  i?(50)  -  0.23 

MB 

R  =  1.26  xi?(50)  -0.44 

R  =  1.50  X  i?(50)  -  0.93 

4dB 

R  =  1.65  xi?(50)  -  1.09 

R  =  2.15  xi?(50)  -2.27 

Table  5.  Regression  result  of  Lambda-Sigma  (1  -  cr)  model  at  each 
X  and  a 

As  illustrated  in  Tables  4  and  5,  it  is  possible  to  estimate 
the  best-fit  R  from  problem  parameters  i?(50),  X,  and  o.  Then 
this  best-fit  R  can  be  used  with  the  Random  Search  formula  below 
to  further  estimate  the  area  search  probability  of  detection 
by  time  t. 

-2RVt 

P(T  <  t)  =  Frit)  =  1  -  e~^. 
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VI .  CONCLUSIONS 


A.  CONCLUSIONS  AND  RECOMMENDATIONS 

There  are  two  primary  contributions  of  this  thesis.  The 
first  is  the  demonstration  that  initial  detection  times  for  area 
search  simulations  using  both  the  Poisson  Scan  and  Lambda-Sigma 
(1-a)  acoustic  detection  models  are  approximately 
exponentially  distributed,  allowing  the  simulation  results  to 
be  closely  approximated  by  the  venerable  Random  Search  formula. 
And  the  second  contribution  is  the  observation  that  the  best-fit 
cookie-cutter  detection  range  used  in  the  Random  Search  formula 
can  be  accurately  predicted  using  the  simulation  model 
parameters  A,  o,  and  the  i?(50)  detection  range. 

In  this  thesis,  it  is  assumed  that  acoustic  signal 
transmission  loss  follows  either  spherical  spreading  or 
cylindrical  spreading.  A  potentially  more  realistic  model  could 
be  developed  by  using  actual  propagation  loss  data. 
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APPENDIX 


A.  MATLAB  CODE 

1.  Poisson  Scan  Model (I) 


Nreps=500;  %number  of  simulation  replications 
tmax=150;  %max.  simulation  time  (hr) 
lx=150;  %search  area  length  in  x  direction  (nm) 
ly=150;  %search  area  length  in  y  direction  (nm) 

V=15;  %searcher  speed  (nm/hr) 

U=5;  %target  speed  (nm/hr) 

lams= (2*V) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
lamt= (2*U) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
laml=2;  %searcher  looking  rate  (1/hr) 
sig=3;  %signal  excess  variance  (dB) 

FOM=60;  %figure  of  merit  (dB) 

R50=10'"  (FOM/15) /1852;  %R(50)  (nm) 
dt=R50/ (2*V) ;  %delta  t  (hours) 

Xs=zeros ( 1 , tmax/dt+1 ) ;  %initialize  x-position  to  zero ( searcher ) 

Ys=Xs;  %initialize  y-position  to  zero ( searcher ) 

Cs=Xs;  %initialize  searcher  course  to  zero 
Xt=Xs;  %initialize  x-position  to  zero (target) 

Yt=Xs;  %initialize  y-position  to  zero (target) 

Ct=Xs;  %initialize  target  course  to  zero 
T=0 : dt : tmax;  %simulation  time  vector 
A=lx*ly;  %search  area 

CumDet=zeros ( 1 ,  tmax/dt+1);  %initialize  cumulative  detection  state 
for  n=l:Nreps  %main  simulation  loop 

xs=rand*lx;  %initial  searcher  and  target  x  and  y  positions 

ys=rand*ly; 

xt=rand*lx; 

yt=rand*ly; 

cs=rand*2*pi;  %initial  searcher  course 
ct=rand*2 *pi ;  %initial  target  course 
t=0;  %set  simulation  time  to  0 
tindex=l;  %initialize  time  index  to  1 
Xs (tindex) =xs;  %save  initial  searcher  x  position 
Ys (tindex) =ys;  %save  initial  searcher  y  position 
Cs (tindex) =cs;  %save  initial  searcher  course 
Xt (tindex) =xt;  %save  initial  target  x  position 
Yt (tindex) =yt;  %save  initial  target  y  position 
Ct (tindex) =ct;  %save  initial  target  course 

Distance=zeros  ( 1 ,  tmax/dt+1);  %initialize  distance  between  target  and 
searcher 

SE=zeros(l,  tmax/dt+1);  %initialize  signal  excess 

Gamma=zeros ( 1 ,  tmax/dt+1);  %initialize  instantaneous  probability  of 
detection 

Detection=zeros ( 1 ,  tmax/dt+1);  %  initialize  detection  vector 
for  t=l:tmax/dt  %inner  loop 
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tindex  =  tindex+1;  %update  simulation  time  index 

if  rand<lams*dt;  cs=rand*2 *pi ;  end 

if  Xs (tindex-1) <0;  cs= (rand-0 . 5) ;  end 

if  Xs (tindex-1 ) >lx;  cs=pi+ (rand-0 . 5) ;  end 

if  Ys (tindex-1) <0;  cs=pi/2+ (rand-0 . 5) ;  end 

if  Ys (tindex-1 ) >ly;  cs=-pi/2+ (rand-0 . 5) ;  end 

if  rand<lamt*dt;  ct=rand*2 *pi ;  end 

if  Xt (tindex-1) <0;  ct= (rand-0 . 5) ;  end 

if  Xt (tindex-1) > (lx) ;  ct=pi+ (rand-0 . 5) ;  end 

if  Yt (tindex-1) <0;  ct=pi/2+ (rand-0 . 5) ;  end 

if  Yt (tindex-1) > (ly) ;  ct=-pi/2+ (rand-0 . 5) ;  end 

Xs (tindex)  =  Xs (tindex-1) +V*dt*cos (cs) ;  %Update  x  and  y  positions 
Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 

Cs (tindex) =cs; 

Xt (tindex)  =  Xt (tindex-1 ) +U*dt*cos (ct) ;  %Update  x  and  y  positions 
Yt (tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct (tindex) =ct; 

Distance (tindex-1 ) =sqrt ( (Xs (tindex-1 ) -Xt (tindex-1 ) )  . ^2+  (Ys (tindex-1 ) -Yt 
(tindex-1 )  )  .''2); 

SE (tindex-1) =FOM-15*loglO (Distance (tindex-1) *1852)  ; 

Gamma (tindex-1) =laml*normcdf (SE (tindex-1) / sig) ; 
if  rand  <=  Gamma (tindex-1 ) *dt; 

Detection (t ; (tmax/dt+1 ) ) =1 ; 

end 

end  %inner  loop  (time  increasing  from  0  to  tmax) 

CumDet  =  CumDet  +  Detection; 
end  %outer  loop  (simulation  replications) 

Probability=CumDet/Nreps ; 

plot (T,  Probability,  'b- ' )  ,  axis ( [0, tmax, 0,1]) 

xlabel ( ' Time ' ,  'Fontsize',  12),  ylabel ( ' CDP ' ,  'FontSize',  12) 
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2. 


Poisson  Scan  Model (II) 


Nreps=500;  %number  of  simulation  replications 
tmax=150;  %max.  simulation  time  (hr) 
lx=150;  %search  area  length  in  x  direction  (nm) 
ly=150;  %search  area  length  in  y  direction  (nm) 

V=15;  %searcher  speed  (nm/hr) 

U=5;  %target  speed  (nm/hr) 

lams= (2*V) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
lamt= (2*U) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
laml=2;  %searcher  looking  rate  (1/hr) 
sig=3;  %signal  excess  variance  (dB) 

FOM=60;  %figure  of  merit  (dB) 

R50=10'' (FOM/15) /1852;  %R(50)  (nm) 
dt=R50/ (2*V) ;  %delta  t  (hours) 

Xs=zeros ( 1 , tmax/dt+1 ) ;  %initialize  x-position  to  zero ( searcher ) 

Ys=Xs;  %initialize  y-position  to  zero ( searcher ) 

Cs=Xs;  %initialize  searcher  course  to  zero 
Xt=Xs;  %initialize  x-position  to  zero (target) 

Yt=Xs;  %initialize  y-position  to  zero (target) 

Ct=Xs;  %initialize  target  course  to  zero 
T=0 : dt : tmax;  %simulation  time  vector 
A=lx*ly;  %search  area 

CumDet=zeros ( 1 ,  tmax/dt+1);  %initialize  cumulative  detection  state 
for  n=l:Nreps  %main  simulation  loop 

xs=rand*lx;  %initial  searcher  and  target  x  and  y  positions 

ys=rand*ly; 

xt=rand*lx; 

yt=rand*ly; 

cs=rand*2*pi;  %initial  searcher  course 
ct=rand*2 *pi ;  %initial  target  course 
t=0;  %set  simulation  time  to  0 
tindex=l;  %initialize  time  index  to  1 
Xs (tindex) =xs;  %save  initial  searcher  x  position 
Ys (tindex) =ys;  %save  initial  searcher  y  position 
Cs (tindex) =cs;  %save  initial  searcher  course 
Xt (tindex) =xt;  %save  initial  target  x  position 
Yt (tindex) =yt;  %save  initial  target  y  position 
Ct (tindex) =ct;  %save  initial  target  course 

Distance=zeros  ( 1 ,  tmax/dt+1);  %initialize  distance  between  target  and 
searcher 

SE=zeros(l,  tmax/dt+1);  %initialize  signal  excess 

Gamma=zeros  ( 1 ,  tmax/dt+1);  %initialize  instantaneous  probability  of 
detection 

Detection=zeros ( 1 ,  tmax/dt+1);  %  initialize  detection  vector 
for  t=l:tmax/dt  %inner  loop 

tindex  =  tindex+1;  %update  simulation  time  index 
if  rand<lams*dt;  cs=rand*2 *pi ;  end 
if  Xs (tindex-1) <0;  cs= (rand-0 . 5) ;  end 
if  Xs (tindex-1 ) >lx;  cs=pi+ (rand-0 . 5) ;  end 
if  Ys (tindex-1) <0;  cs=pi/2+ (rand-0 . 5) ;  end 
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if  Ys (tindex-1 ) >ly;  cs=-pi/2+ (rand-0 . 5) ;  end 
if  rand<lamt*dt;  ct=rand*2 *pi ;  end 
if  Xt (tindex-1) <0;  ct= (rand-0 . 5) ;  end 
if  Xt (tindex-1) > (lx) ;  ct=pi+ (rand-0 . 5) ;  end 
if  Yt (tindex-1) <0;  ct=pi/2+ (rand-0 . 5) ;  end 
if  Yt (tindex-1) > (ly) ;  ct=-pi/2+ (rand-0 . 5) ;  end 

Xs (tindex)  =  Xs (tindex-1) +V*dt*cos  (cs) ;  %Update  x  and  y  positions 
Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 

Cs (tindex) =cs; 

Xt (tindex)  =  Xt (tindex-1 ) +U*dt*cos (ct) ;  %Update  x  and  y  positions 
Yt (tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct (tindex) =ct; 
if  rand<laml*dt 

Distance  (tindex-1 )  =sqrt  (  (Xs  (tindex-1 )  -Xt  (tindex-1 )  )  .  '^2  +  (Ys  (tindex-1 )  -  Yt 
(tindex-1 ) )  .  ^2 )  ; 

SE (tindex-1) =F0M-20*logl0 (Distance (tindex-1) *1852)  ; 

Gamma (tindex-1) =normcdf (SE (tindex-1) / sig)  ; 
if  rand  <=  Gamma (tindex-1 ) ; 

Detection (t : (tmax/dt+1 ) ) =1 ; 

end 

end 

end  %inner  loop  (time  increasing  from  0  to  tmax) 

CumDet  =  CumDet  +  Detection; 
end  %outer  loop  (simulation  replications) 

Probability=CumDet/Nreps ; 

plot (T,  Probability,  'b- ' )  ,  axis ( [0, tmax, 0,1]) 

xlabel ( ' Time ' ,  'Fontsize',  12),  ylabel ( ' CDP ' ,  'FontSize',  12) 


42 


3. 


Lambda-Sigma  (A  —  a)  Model 


Nreps=500;  %number  of  simulation  replications 
tmax=150;  %max.  simulation  time  (hr) 
lx=150;  %search  area  length  in  x  direction  (nm) 
ly=150;  %search  area  length  in  y  direction  (nm) 

V=15;  %searcher  speed  (nm/hr) 

U=5;  %target  speed  (nm/hr) 

lams= (2*V) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
lamt= (2*U) /sqrt (lx*ly) ;  %searcher  course  change  rate  (1/hr) 
lamda=2;  %duration  of  each  step  (1/hr) 
sig=3;  %height  of  each  step  (dB) 

FOM=60;  %  figure  of  merit  (dB) 

R50=10'' (FOM/15) /1852;  %R(50)  (nm) 
dt=R50/ (2*V) ;  %delta  t  (hours) 

Xs=zeros ( 1 , tmax/dt+1 ) ;  %initialize  x-position  to  zero ( searcher ) 
Ys=Xs;  %initialize  y-position  to  zero ( searcher ) 

Cs=Xs;  %initialize  searcher  course  to  zero 
Xt=Xs;  %initialize  x-position  to  zero (target) 

Yt=Xs;  %initialize  y-position  to  zero (target) 

Ct=Xs;  %initialize  target  course  to  zero 
SP=Xs;  %initialize  step  function 
T=0 : dt : tmax;  %simulation  time  vector 
A=lx*ly;  %search  area 

CumDet=zeros ( 1 ,  tmax/dt+1);  %initialize  cumulative  detection  state 
for  n=l:Nreps  %main  simulation  loop 

xs=rand*lx;  %initial  searcher  and  target  x  and  y  positions 

ys=rand*ly; 

xt=rand*lx; 

yt=rand*ly; 

cs=rand*2 *pi ;  %initial  searcher  course 
ct=rand*2*pi;  %initial  target  course 
sp=randn*sig;  %initial  step 
t=0;  %set  simulation  time  to  0 
tindex=l;  %initialize  time  index  to  1 

Xs (tindex) =xs;  %save  initial  searcher  x  position 

5save  initial  searcher  y  position 

jsave  initial  searcher  course 

jsave  initial  target  x  position 

5save  initial  target  y  position 

jsave  initial  target  course 

5save  initial  step 


Ys (tindex) =ys; 
Cs (tindex) =cs; 
Xt (tindex) =xt; 
Yt (tindex) =yt; 
Ct (tindex) =ct; 


SP (tindex) =sp; 

Distance=zeros  ( 1 ,  tmax/dt+1) 
searcher 

SE=zeros(l,  tmax/dt+1);  %initialize  signal  excess 


%initialize  distance  between  target  and 


Detection=zeros ( 1 ,  tmax/dt+1);  %initialize  detection  vector 
for  t=l:tmax/dt  %inner  loop 

tindex  =  tindex+1;  %update  simulation  time  index 
SP (tindex) =SP ( tindex- 1 ) ; 

if  rand<lamda*dt;  SP (tindex) =randn*sig;  end 
if  rand<lams*dt;  cs=rand*2 *pi ;  end 
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if  Xs (tindex-1) <0;  cs= (rand-0 . 5) ;  end 
if  Xs (tindex-1 ) >lx;  cs=pi+ (rand-0 . 5) ;  end 
if  Ys (tindex-1) <0;  cs=pi/2+ (rand-0 . 5) ;  end 
if  Ys (tindex-1 ) >ly;  cs=-pi/2+ (rand-0 . 5) ;  end 
if  rand<lamt*dt;  ct=rand*2 *pi ;  end 
if  Xt (tindex-1) <0;  ct= (rand-0 . 5) ;  end 
if  Xt (tindex-1) > (lx) ;  ct=pi+ (rand-0 . 5) ;  end 
if  Yt (tindex-1) <0;  ct=pi/2+ (rand-0 . 5) ;  end 
if  Yt (tindex-1) > (ly) ;  ct=-pi/2+ (rand-0 . 5) ;  end 

Xs (tindex)  =  Xs (tindex-1) +V*dt*cos (cs) ;  %Update  x  and  y  positions 
Ys (tindex)  =  Ys (tindex-1) +V*dt*sin (cs) ; 

Cs (tindex) =cs; 

Xt (tindex)  =  Xt (tindex-1 ) +U*dt*cos (ct) ;  %Update  x  and  y  positions 
Yt (tindex)  =  Yt (tindex-1 ) +U*dt*sin (ct) ; 

Ct (tindex) =ct; 

Distance  (tindex-1 )  =sqrt  (  (Xs  (tindex-1 )  -Xt  (tindex-1 )  )  .  "'2  +  (Ys  (tindex-1 )  -  Yt 
(tindex-1 )  )  .''2); 

SE (tindex-1) =FOM-15*loglO (Distance (tindex-1) *1852) ; 
if  SP (tindex-1)  <=  SE (tindex-1 ) ; 

Detection (t :  (tmax/dt+1 ) ) =1 ; 

end 

end  %inner  loop  (time  increasing  from  0  to  tmax) 

CumDet  =  CumDet  +  Detection; 
end  %outer  loop  (simulation  replications) 

Probability=CumDet/Nreps ; 

plot (T,  Probability,  'b- ' )  ,  axis ( [0, tmax, 0,1]) 

xlabel ( ' Time ' ,  'Fontsize',  12),  ylabel ( ' CDP ' ,  'FontSize',  12) 
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